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Recent advances in our understanding of the higher- wave quarkonia have generated much interests in quarko- 
nium physics. However most are devoted to inclusive decays and productions. Experimental data of several 
two-body exclusive decay channels of P-wave charmonia such as irn and pp are available and some have recently 
been re-measured by the BES collaboration. It is not clear from the outset that color octet is needed for these 
exclusive channels. Indeed only color singlet has been used in the past and reasonable agreement with data was 
found. Contrary to these old results, we provide theoretical arguments for the inclusion of color octet and perform 
explicit calculations to back this up. 



1. Introduction 

Quarkonia are special hadronic systems in that 
they have the mass as an intrinsic large scale. 
This permits factorization, the separation of the 
hard from the soft scale physics. Therefore once 
supplemented with non-perturbative matrix ele- 
ments in inclusive processes or hadronic wave- 
functions in exclusive processes, physical pro- 
cesses involving heavy quarkonia can be calcu- 
lated perturbatively. In this talk, we would like 
to re-examine some exclusive charmonium de- 
cay channels. There are several reasons for do- 
ing this. First being non-relativistic, the heavy 
quark-antiquark system is simpler than most of 
the lighter hadrons and therefore is good for the 
study of color confinement in this special case. 
If sufficiently understood, it provides a testing 
ground for our understanding of wavefunctions of 
selected lighter hadrons. Then recently several 
decay channels of P-wave charmonia have been 
measured or re-measured at the Beijing Spec- 
trometer by the BES collaboration so there are 
new data available. Third if one looks up the ex- 
isting calculations, it is immediately evident that 
they are rather out-of-date. Not only the QCD 
parameters have changed, the hadronic wavefunc- 
tions used have also been updated and better 
understood. Last but not least there is also 
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the emergence of the color octet higher states 
in quarkonium, their significance was revealed in 
|j],|2J . The simplest decays are the two-body chan- 
nels and so we will look at these of the P-wave \ J 
system. 

2. Two-body P-wave charmonium decay 

A well-known scheme for calculating exclusive 
processes is the hard scattering approach (sHSA) 
of Brodsky and Lepage || which became a kind 
of standard so that is what the small "s" stands 
for in the acronym. This scheme relies on the 
presence of a large momentum transfer to vali- 
date factorization so that the probability ampli- 
tude M of the process in question can be written 
as a convolution of the hadronic distribution am- 
plitudes 4>h of the hadrons involved and a hard 
perturbative part Tu which serves as the conduit 
of the large momentum flow in the Feynman di- 
agram between the initial and final hadrons. For 
a two-body decay process resulting in hadron h\ 
and /12, this can be summarized in the following 
form 



M 



(1) 



fxAxJ ® T H ® hAhi ® fh 2 4>h 2 

The fh's are the decay constants of the respec- 
tive hadrons. The convolution here is over the 
light-cone momentum fractions. The partial de- 
cay width in terms of the above amplitude and 
the heavy charmonium mass M is essentially 
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Figure 1. The graphs for calculating the color 
singlet contribution for the \J decay into tttt. 
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Figure 2. The color singlet graphs for the decay 
into pp. 



For definiteness we will examine the \J decay into 
tttt and pp. Both of these have been measured and 
data are available from both Particle Data Group 
g] and BES @. 

The calculations of the decay into the above 
two channels are quite reasonable with only two 
diagrams in the former and four in the latter pro- 
cess and so are not too complicated [BUT]. These 
are shown in Fig. [I] and Fig. g. Using the sHSA 
and some rather well-known quantities such as 
the /„•, \J radial wavefunction i?p(0), m c etc. 
the partial widths can be obtained. One gets the 
results shown in Table w. Comparing the partial 
widths T^ 1 ' coming from the cc state to those from 
experiments, one can see that the former are well 
below the data. 

In order to make sure that these small val- 
ues are not the result of the imperfection of the 
scheme used. They can be calculated again us- 
ing an improved scheme (mHSA) of Sterman et 
al [|8|,|9|. This scheme comes from modifying the 
sHSA by introducing transverse momentum ef- 
fects, radiative correction in the form of Sudakov 
factor S and utilizing whole hadronic wavefunc- 
tions iph- The probability amplitude can now be 



Table 1 

Partial decay widths obtained within the sHSA 

scheme. 
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Table 2 

Partial decay widths obtained within the mHSA 

scheme. 
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summarized as 

M ~ tp XJ <g) T H (a s ) <8> il> hl <8> Vh 2 <8> expl-S*} .(3) 

The convolution is now over not just light-cone 
fractions but also internal transverse momenta. 
Furthermore the coupling a s is part of the convo- 
lution whose scale is determined dynamically by 
the momentum flow, whereas in sHSA it is merely 
a constant. The importance of this was discussed 
in [0,|l(| . The results shown in Table || agree es- 
sentially with those in Table pi on the fact that 
the theoretical results are too small. Of course, 
one could exploit the uncertainties in the various 
parameters to try to enlarge the calculated results 
but this proved to be fruitless. They remain be- 
low the data by factors of two or more pMO,[12| . 



3. What is the source of the problem? 

The previous section showed that it was not 
possible to explain the experimental data from 



the decay of the P-wave cc system. In the search 
for the source of the problem, one cannot help 
but wonder if color octet is the answer. After all 
color octet seems these days to be inextricably 
intertwined with quarkonium physics in inclusive 
processes ever since it was first introduced in |jl],|) 
to cancel the infrared divergence found in the in- 
clusive xj decay ]1^] . Strangely in the world of 
exclusive quarkonium processes, color octet does 
not seem to exist at all. In fact it has been totally 
neglected. A simple argument easily reveals that 
color octet must also be important in exclusive 
decays, albeit it is not sufficient to claim that it 
is so in the two-body channels. The argument re- 
lies on the fact that the sum of the partial widths 
of all channels must be equal to the total inclu- 
sive hadronic width which can be more precisely 
expressed in the form 



r — r 



/ -L i ^ incl. 

i£ channels 
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If none of the I\ requires the color octet contri- 
bution, then there would be a contradiction be- 
cause it is known to be of importance and needed 
on the right-hand-side (r.h.s). However the con- 
tradiction can be avoided provided color octet is 
needed in at least one channel, it need not be the 
two-body channels that we are interested in. So a 
priori for our problem at hand, it is not clear that 
color octet is the answer. Furthermore there is 
the usual folklore that says higher state contribu- 
tions are suppressed in large momentum transfer 
processes which seems to oppose the inclusion of 
the octet contribution. 

4. The large mass dependence of the decay 
amplitudes 

The pros and cons of color octet as the answer 
to our problem discussed in the last section makes 
the issue very unclear. To come to a resolution, 
it is necessary to look closely at the charmonium 
system. The hadronic decay is via annihilation 
into gluons. The cc pair has to come close to- 
gether for this at a small distance I ~ l/M be- 
cause of the large charmonium mass M. For a 
P-wave charmonium, the L = 1 orbital angular 
momentum tends to force the pair apart there- 



fore it is harder for the pair to annihilate. This is 
manifested in the vanishing of the P-wave wave- 
function at the origin tpp(0) — 0. Comparing a 
S-wave to a P-wave system, assuming all else be- 
ing equal, the only difference in magnitude of the 
decay amplitude lies entirely with the wavefunc- 
tions at small distance. Because the P-wave wave- 
function vanishes there, its derivative enters the 
decay amplitude in its place. The following shows 
the relevant quantities that enter the amplitudes 
and the same after transforming into momentum 
space. 



S : Ml ~ 0) 

P: ip P {l ~ 0) ~ lip' P (l ~ 0) 






This shows the P-wave quantity carried with it a 
power of l/M so a P-wave charmonium is sup- 
pressed at the valence level already by M in com- 
parison to a S-wave. In view of this power sup- 
pression, it is perhaps useful to examine and com- 
pare the M dependence of the decay amplitudes. 
In the calculations of exclusive process, one has 
to deal all too frequently with decay constants of 
the hadrons involved. So it should be more rele- 
vant to look at the decay amplitudes in terms of 
fh's and M. The sHSA scheme is therefore more 
convenient for this purpose. Since the strongest 
dependence on M is powerlike, it is useful to look 
at dimensions. From Eq. (0) M. is of mass di- 
mension one so we have the following equation of 
dimensions 



[M] = [mass] 1 = [f XJ ][f hl ][h 2 )[T H ] . 



(5) 



The r.h.s. shows all the quantities that carried a 
mass dimension. The hard part Tu has some hid- 
den power dependence of the form M~ p for some 
number p. Together with the /h's, they make up 
for the dimension of M.. The M dependence of 
M will be known, once the mass dimensions of 
the fhS are determined. Using the relation of ft 
to the associated wavefunction iph 
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and the normalization of iph 
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one can deduce the dimension of fh- It can be 
summarized by 



[i'h] = [mass 



l-N 



[f h ] = [mass] N - 1+L (8) 



where N is the number of constituents in the 
hadron h described by iph and L is the orbital 
angular momentum of iph- The decay constants 
for our hadrons are then 

[U] = [mass] 1 , [f p ] = [/«] = [/£>] = [mass] 2 

For the XJ decays, the color singlet and octet 
amplitude therefore have the following M depen- 
dence 
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For both decay channels, singlet and octet con- 
tribution have the same power dependence in M. 
There is no large M suppression of the octet in 
comparison to the singlet. The usual suppression 
of the higher state is nullified by the suppression 
of the P-wave valence state due to angular mo- 
mentum. 

To elucidate this further, we can turn to the 
same decay channels of the S-wave J/ip. Now the 
J/ip decay constants have the dimensions 
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The octet constant has a dimension of three in- 
stead of two is because the three-body wavefunc- 
tion must have L — 1 to make up the right quan- 
tum numbers for J/ip as explained in W]. The 



J/ip decay amplitudes have the dependence 
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So the octet contribution is suppressed by M~ 2 
in both cases in the amplitude. One can therefore 
legitimately neglect the color octet contributions 
in J/ip but not in \J- 

5. Color octet contributions 

In this section, we briefly outline and discuss 
the calculation of the octet contributions. The 
calculation can be done in either sHSA or mHSA 
scheme in principle, in practice it is easier to use 
sHSA than mHSA but the latter has important 
advantages discussed in [MITM!^ J14J . To some de- 
gree the choice of scheme depends also on the final 
hadrons. The simpler rnr is manageable in mHSA 
but not the pp channel which is too complicated. 
It is better to work in sHSA for calculating the 
latter. 

Since the color octet wavefunctions are not 
known, they have to be constructed. To do this 
the associated decay constants f^J must be de- 
termined by fitting and the light-cone fractions' 
distribution fixed in a sensible manner. Numer- 
ically the latter has to be checked to see if it is 
sensible. Then there is the presence of the con- 
stituent gluon from the octet state of the char- 
monium which has to be considered. The usual 
valence color singlet cc state is straight forward 
to deal with since they cannot appear in the final 
state and easily eliminated via annihilation. The 
constituent gluon on the other hand can be a part 
of the final hadron by becoming a constituent of a 
higher state or it can end in the hard perturbative 
part Th- It was eventually decided that the dom- 
inant contribution would be for it to end on Th- 
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Figure 3. The groups of graphs that contribute 
to the decay into 7T7T. 



The contribution involving the higher state of a 
final hadron would be suppressed both by large 
momentum flow within HSA and the presence of 
the additional higher state. To complete the cal- 
culation of Tjj , all possible feynman graphs must 
be found. Unlike the color singlet contribution, 
there are now a lot more diagrams. Because of C- 
parity, the color singlet state must annihilate into 
two or more gluons. The cc in the color octet on 
the other hand is protected from the C-parity con- 
straint by the presence of the constituent gluon 
and can therefore annihilate through one gluon. 
This leads to quite a number of diagrams. Addi- 
tionally by allowing the constituent gluon to end 
in the hard part, it must be attached to all possi- 
ble allowed positions of Th ■ These all contribute 
to the number of graphs. In the end, they can 
be organized into groups and calculated by com- 
puter. Some examples are shown in Fig. H and 
Fig. [|. The first figure here is from M . The con- 
stituent gluon is yet to be included in the latter. 
The graphs shown in this one form the basis from 
which the groups are to be generated. It remains 
to convolute the (j)hS in the sHSA or iph in the 
mHSA and Th together to get the partial widths 
for the different two-body channels. The details 
for the calculation of the \J decay into tttt are to 
be found in B,Q and those for \J decay into pp 
are in |]|l2| . 
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Figure 4. The groups that contribute to the decay 
into pp. The constituent gluon is not drawn here. 



6. Combining color singlet and octet con- 
tributions 

The arguments in Sec. Hshow that for P-wave 
charmonium, color octet and singlet are of equal 
weight and must both be included. With the cal- 
culation of the singlet contribution in Sec. |2| and 
the octet contribution in Sec. g, the two can now 
be combined to give, from a theoretical point of 
view, the proper partial widths of xj- The results 
are tabulated in Table ||Q. Notice that the exper- 
imental measurements are different between the 
PDG and BES data. It is up to the experimental- 
ists to reconcile this in the future. The important 
point here is that by including the octet contribu- 
tions, the theoretical widths are in much better 
shape for explaining the measured data and also 
in accordance with theoretical expectation. With 
hindsight it is not entirely surprising that color 
octet is required in the two-body decay channels 
since in inclusive decays it is needed already at 
the leading order for xi an d next-to-leading or- 
der for xo and %2- 



2 These results for the decay into proton-antiproton are 
final. The y sup ersede all previously reported preliminary 
values in 
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Table 3 

Singlet and singlet-octet combined partial widths. 
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105.0 ± 30 64.0 ± 21.0 
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56.27 
154.19 


75.68 ± 10.5 37.84 
200.00 ± 20.0 118.00 


7. Summary 





In this talk, we showed with new calculations 
that valence color singlet state of P-wave char- 
monia alone is insufficient to account for experi- 
mental measurements of the partial decay widths 
in contradiction to those done in the past pre- 
dating the emergence of the color octet state 
in quarkonium physics. Normally higher states 
should be suppressed but this was shown to be 
circumvented because of angular momentum sup- 
pression of the valence state with L = 1 which 
brought it down to the same level of the octet 
state. Power dependence in the heavy charmo- 
nium mass arguments revealed that it was in- 
consistent to perform P-wave charmonium calcu- 
lations without including the color octet. This 
was eventually supported by explicit calculations. 
The anormaly of the situation in charmonium 
physics between inclusive and exclusive processes 
pointed out was thus reconciled. The arguments 
given here can be straight forwardly generalized 
to even higher-wave quarkonia. Calculations in- 
volving the latter require the inclusion of not only 



the next higher states but also those above them. 
Therefore one has to be very careful when dealing 
with quarkonia whose valence states have non- 
zero orbital angular momentum. 
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